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Abstract —A counter braid (CB) is a novel counter architecture 
introduced by Lu et al. in 2007 for per-flow measurements on 
high-speed links. CBs achieve an asymptotic compression rate 
(under optimal decoding) that matches the entropy lower bound 
of the flow size distribution. Spatially-coupled CBs (SC-CBs) 
have recently been proposed. In this work, we further analyze 
single-layer CBs and SC-CBs using an equivalent bipartite graph 
representation of CBs. On this equivalent representation, we show 
that the potential and area thresholds are equal. We also show 
that the area under the extended belief propagation (BP) extrinsic 
information transfer curve (defined for the equivalent graph), 
computed for the expected residual CB graph when a peeling 
decoder equivalent to the BP decoder stops, is equal to zero 
precisely at the area threshold. This, combined with simulations 
and an asymptotic analysis of the Maxwell decoder, leads to the 
conjecture that the area threshold is in fact equal to the Maxwell 
decoding threshold and hence a lower bound on the maximum 
a posteriori (MAP) decoding threshold. Finally, we present some 
numerical results and give some insight into the apparent gap of 
the BP decoding threshold of SC-CBs to the conjectured lower 
bound on the MAP decoding threshold. 

I. Introduction 

Recently, Lu et al. proposed a novel counter architecture, 
inspired by sparse graph codes, for measuring network flow 
sizes, nicknamed Counter Braid (CB) ID-0- CBs use less 
memory space than other approximate counting techniques, 
since the flow sizes are compressed on-the-fly. They are 
indeed asymptotically optimal (under some mild conditions) 
0, i-e., the average number of bits needed to store the 
size of a flow tends to the information-theoretic limit (under 
maximum-likelihood decoding) when the number of flows 
goes to infinity. Furthermore, they are characterized by a 
layered structure which can be described by a graph. This 
allows for low-complexity decoding using a message passing 
(belief propagation (BP)) decoding algorithm. In general, good 
performance can be achieved with a small number of layers. 

Spatially-coupled CBs (SC-CBs) were recently introduced 
by Rosnes in 0, with improved BP decoding thresholds as 
compared to uncoupled CBs. It was numerically shown in 
0 that the BP decoding threshold of single-layer SC-CBs 
converges to a fixed value. 

In this paper, we further analyze single-layer CBs and SC- 
CBs. We consider an equivalent bipartite graph representation 
of CBs, with identical finite-length performance and asymp¬ 
totic behavior to that of CBs decoded on the original bipartite 
graph. On this equivalent representation, we show that the 
potential threshold, introduced by Yedla et al. in 0, and the 
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Fig. 1. The graph of a two-layer CB. (a) the bipartite graphs Q\ and 
( b ) the combined graph. 

area threshold are equal. We also show that the area under 
the extended BP (EBP) extrinsic information transfer (EXIT) 
curve (defined for the equivalent graph), computed for the 
expected residual CB graph when a peeling decoder equivalent 
to the BP decoder stops, is equal to zero precisely at the area 
threshold. This, combined with simulations and an asymptotic 
analysis of the Maxwell decoder, leads to the conjecture that 
the area threshold is equal to the Maxwell decoding threshold 
and hence a lower bound on the maximum a posteriori (MAP) 
decoding threshold. 

Interestingly, when coupling the original (or the equivalent) 
graph, there is a remaining gap between the potential threshold 
and the BP decoding threshold of SC-CBs even in the limit 
of large coupling chain length and smoothing parameter. This 
gap seems to be fundamental and due to the fact that the flow 
node update rule for even and odd iterations is different. 

In this text, all proofs are omitted due to lack of space. An 
extended version containing proofs will be published on arXiv. 

II. Counter Braids 

A CB is a counter architecture consisting of L > 1 layers. 
At layer l £ {1,..., L}, it has mi counters of depth di bits, 
with Wj < rrij for i > j. The number of distinct flows to be 
counted is denoted by too. 

The Z-th layer of a CB can be represented by a bipartite 
graph Qi = Qi(Ti U Ci,£i), where J) (of size to/_i) denotes 
the set of flow nodes, and Ci (of size mi) denotes the set of 
counter nodes. The set of edges of the graph is denoted by 
Si. The Z-th layer is connected to the (Z — l)-th layer by a 
bijective mapping E)(/) = c on the set of flow nodes of the 
Z-th layer, where / £ Ti and c £ Cz_i, Z = 2. ... . L. For 
notational convenience, the neighborhood of a node a (either 
a counter node c or a flow node /) is denoted by 6(a). 

When a flow is encountered (for instance, on a high¬ 
speed link), all connected counter nodes of the first layer 















are incremented modulo 2 dl . If a counter node c of the first 
layer overflows, all connected counter nodes of the second 
layer (formally the counter nodes in the set J(E^' 1 (c))), are 
incremented modulo 2 d2 . Furthermore, if a counter node in 
the second layer overflows, all connected counter nodes of 
the third layer are incremented modulo 2“ 3 . This process is 
repeated for each level until we reach the final layer L. We 
denote by </>(c) the final value of a counter node c prior to 
decoding, and by </>(/) the estimated value (after decoding) of 
a flow node /. 

An example of a two-layer CB is shown in Fig. [Q where 
flow nodes are represented by empty circles and counter 
nodes by filled squares. Fig. Uf a) shows the bipartite graphs 
Qi and Q 2 of the two-layer CB, while Fig. [TJ fo) depicts an 
equivalent graph where a flow node of the second layer and 
its corresponding counter node of the first layer in Fig. Oja) 
are represented using a single combined counter node due to 
the bijective mapping £ 2 (•)■ 

In the following, we will assume that the bipartite graphs 
Qi, l = 1are left-regular, i.e., \S(f)\ = ki, for some 
integer ki > 2, and for all / g J;. The assignment of graph 
sockets of flow nodes to graph sockets of counter nodes, i.e., 
the connections in the graph, is done in a random fashion. This 
means that, asymptotically as the number of flow nodes tends 
to infinity (while is kept fixed), the distribution of the 
fraction of edges connected to a given counter node approaches 
a Poisson distribution (TJ. We can now formally define the flow 
and counter node degree distributions. For counter nodes, we 
assume the asymptotic Poisson distribution. Let Li(z ) = z kl , 
ki > 2 , and 


R l {z) = Y,R? zi 

i =0 


2 = 0 
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denote the node-perspective flow node and counter node 
degree distributions of the Z-th layer, respectively. Here, 7 1 = 
m ‘ m 1 [ k ‘ is the average counter node degree. The corresponding 
edge-perspective degree distributions are 
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respectively. Finally, we denote by /3i = the fraction of 
the number of counter and flow nodes of the l -th layer. 


A. Belief Propagation Decoding 

CBs can be decoded using a BP decoding algorithm m 
on the corresponding graph, starting with the right-most layer. 
After decoding the Z-th layer, represented by the bipartite graph 
Qi, the counter nodes of the (/ — l)-th layer are updated based 
on the mapping £/(•). 

Consider the decoding of the /-th layer in more detail. We 
define the messages exchanged in the BP decoding algorithm 
as follows. Let M/4-c ^ where R is the field of real numbers 


and (/, c) £ £1 , denote the message sent from flow node / to 
counter node c for some Z, during the £-th iteration of the 
BP decoding algorithm. Likewise, let € R, (/, c) £ £ 1 , 
denote the message sent from counter node c to flow node / 
during the £-th iteration of the algorithm. The counter node 
and flow node update rules (for £ = 1 ,..., £ max ) are given as 
follows m, 


fpcX f = max { <t >( C ) - H 

fes(c)\f 
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where / m i n is the minimum flow size (from the flow size 
distribution) of a flow node, £ max is the maximum number 
of iterations, and pf\ c = /min for all (/, c) £ £ 1 . The final 
estimate of the flow sizes of the flow nodes is according to 


fcf) 


min cg5(/) if 4 iax is odd 

U J ) 

max ce , 5 ( / ) if tmax is even 


After decoding the Z-th layer (starting with the L- th layer), 
the counter nodes of the (l — l)-th layer are updated as follows 
(for all / £ Ti), 


mu)) ^ kf) -2 di -' + mu))- 


Decoding proceeds layer-by-layer until the first layer is de¬ 
coded. 

In the rest of the paper we assume a single-layer system 
and frequently omit, for notational convenience, the sub¬ 
script/superscript l. Also, we assume that the depth of the 
counters d\ is large enough so that they do not overflow. 


III. Asymptotic Analysis of Single-Layer Counter 
Braids 


Consider a single-layer CB with m 0 flow nodes and m 1 
counter nodes of depth d±. When m 0 —> 00, the asymptotic 
performance of CBs can be analyzed by means of density 
evolution (DE). Denote by and y <f: i the error probability 
of an outgoing message from a flow node and a counter node, 
respectively, at the f-th iteration. The flow node and counter 
node DE updates at iteration £ are described by the equations 


x m = f(v m f, = S (*<'-») , 

where 

~ \y k ~ l , if £ is odd 

f{V\z) = \ k—l -to- 

I e • y , if £ is even 

<?(x) = 1 -p(l -x) 


(3) 
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and e is the probability of observing a flow of size strictly 
larger than / min . 

In the following, we denote by X = [0,1] and y = [0,1] 
the set of possible values for x and y, respectively, and by £ 
the set of possible values for e. 








For a given flow size distribution, or more precisely, for a 
given e, a relevant parameter for the performance of single¬ 
layer CBs is the number of counters per flow, /? = = 

The average number of bits needed to represent a flow is 
therefore /3di, hence 3 is directly related to the compression 
rate. In particular, we are interested in the minimum value of 
f) so that decoding is successful, 

£bp = 0 B p(e) = inf {/3 > 0 | ar«°°>(/3, e) = o} . 


1) The first derivative of f(y; e) exists and is continuous on 
y x £, the first derivative of g(x) exists and is continuous 
on X, 

2 ) f(y ; e) is nondecreasing in both y and e, 

3) g{x) is strictly increasing in x, and 

4) the second derivative of g(x) exists and is continuous 
on X. 

It is easy to show that the DE updates in ©-(dll satisfy 
these properties. 


Alternatively, we can analyze the asymptotic behavior of CBs 
by fixing f3 and finding the maximum value of e such that 
decoding is successful, 

£bp = £bp(/3) = sup je G £ | a; (oo) (/?,e) = oj 

since for a fixed k and any e it follows that £Bp(/3Bp(e)) = 
/3gp(/?Bp(e)) = £■ In this case, £bp has a similar meaning 
as the BP decoding threshold for low-density parity-check 
(LDPC) codes over the binary erasure channel, where e can 
now be interpreted as the channel parameter. 

Combining two successive iterations into a single one, we 
can write the DE in an equivalent form as 

yM =g(x^-V), (6) 

where 

f(y,e) = e-y k ~ 1 (7) 

9{x) = 1 - P (l - (1 - P (1 - a;)) fe_1 ) • (8) 

Clearly, the DE recursions in ©-© and ©-© give the 
same BP decoding threshold. 

There is a way to construct an equivalent graph representa¬ 
tion of single-layer CBs (and a corresponding message passing 
decoding algorithm) which will give identical finite-length 
performance to that of the original graph, and asymptotically 
the DE recursion in ©-©. The new graph has the same 
flow nodes as the original graph and the same number of 
counter nodes, but each counter node represents now a tree 
of depth three grown from the corresponding counter node in 
the original graph. The overall construction is omitted due to 
lack of space. Thus, in the following, we can equivalently base 
our analysis on the recursion in ©-©. 

Lemma 1: The functions f(y; e) (for a fixed value of e >0) 
and g(x), k > 2 , are strictly increasing in y and x, respectively. 

As a result of Lemma [T| it follows that the DE recursion in 
©-© (and thus the DE recursion in ©-©) for uncoupled 
CBs converges to a fixed-point. The fixed-point DE equation 
for x = x(°°\p, e) is 

x = f{g{x)). (9) 

Decoding is successful if the fixed-point is x^°°^ (/3, e) = 0. 

We will need the definition of an admissible system 0 . 

Definition 1 (Admissible System): An admissible system is 
a system where the functions f(y,e) and g(x) satisfy the 
following properties. 


A. Extended Belief Propagation Extrinsic Information Trans¬ 
fer Curve 

The EBP EXIT curve of a single-layer CB is given in 
parametric form by 

(e, h BBP ) = (e(x), (1 - p( 1 - (1 - p{ 1 - x)) k ~'))% ( 10 ) 
where 


W (l-p(l-(l-p(l-aO)*-i))*-i 

is the solution of © for e, and x C X. The curve is a trace 
of all fixed-points of the DE recursion in © (both stable and 
unstable). 

Definition 2 (Area Threshold): Let (e(a;*), ft EBP (a:*)) be a 
point on the EBP EXIT curve ft, EBP of a single-layer CB such 

h EBP (x) de(or) 

and there exist no x' € (x* , 1 ] such that e{x') = e{x*). Then, 
the area threshold, denoted by e, is defined as 6 = e(tr*). 

B. Potential Function, Potential Threshold, and Area Thresh¬ 
old 

Since the DE recursion in ©-© describes an admissible 
system (see Definition ©. we can define a corresponding 
potential function as in 0 for the uncoupled system. 

Definition 3: The (uncoupled) potential function U(x\ e) of 
the system defined by the functions / and g from © and ©, 
respectively, is given by 

U{x;e) = xg{x) - f g(z) dz - f f(z;e) dz 
Jo Jo 

= xg( x )-J^ g(z) dz - ( -g{x) k . 

Following © Def. 32], we make the following definitions, 

’L(e) = m.\nU{x\e), 

x$LX 

X*(e) = {x G X | U(x\ e) = ’L(e)}, and 
x*(e) = maxX*(e). 

Now, we can define the potential threshold e* as 0 Def. 35] 

e* = sup{e e £ \ x*(c) = 0}. (11) 

The following theorem shows that the area threshold and 
the potential threshold are equal, e = e*. 

Theorem 1: The area threshold from Definition [2] is equal 
to the potential threshold from (ITTl i. 




IV. Maxwell Decoder 


V. Spatially-Coupled Counter Braids 


A Maxwell decoder can be constructed for CBs as for 
LDPC codes, and it can be analyzed asymptotically using 
DE on the equivalent graph representation mentioned in Sec¬ 
tion Hn] An important point to consider is that the Maxwell 
decoder is in general not a MAP decoder as for LDPC codes, 
since the flow size distribution is in general nonuniform. Thus, 
the Maxwell decoding threshold (on e) is in general a lower 
bound on the MAP decoding threshold. 

Now, note that the BP decoder with the update rules in ([TJ- 
© will always stop , i.e., at some point the flow size estimates 
of the flow nodes will be the same at iterations £ and £ — 2 for 
some £. When the decoder stops, either the flow size estimates 
for a flow node at iterations l and £ — 1 are the same (i.e., we 
have convergence (upper and lower bounds are the same and 
thus we have the correct value for the flow size)), or we have 
oscillation (i.e., the estimates for iterations £ and i — 2 are the 
same, but the estimates for iterations £ and £ — 1 are not the 
same). 

The BP decoder can be turned into a peeling decoder in 
the following way. Run the BP decoder and in addition apply 
the following two peeling rules: 1) In case |<5(c) ] = 1 for a 
counter node c, then remove c and the connected flow node 
/ £ (5(c) and all its edges from the graph. Decrease the values 
of the counter nodes of 5(f) \ c by the value of c. 2) For 
odd iterations, if a message from a counter node c to a flow 
node / is equal to / m j n , then remove the flow node / and all 
its edges from the graph. Decrease the values of the counter 
nodes of 5(f) by / m ; n . The BP and the peeling decoders are 
equivalent, in the sense that the set of converged flow nodes 
for the BP decoder is exactly equal to the set of peeled flow 
nodes for the peeling decoder. 

Theorem 2: The EBP EXIT curve (defined for the equiva¬ 
lent graph as in (ITOt ) for the expected residual CB graph when 
the peeling decoder stops is given in parametric form by 

(e, h EBF ) = (■£(z ; x), (1 - p( 1 - z; x)) k ), 


where 

e(z-,x) 


z 

(1 - p( 1 - Z-tX))*- 1 


and p(z; x) = 1 — 


g(x ~ zx) 
g(x) 


and where x = x(e) is the largest fixed-point of the DE 
recursion in ©, for a given e. Furthermore, for e = e (the area 
threshold from Definition^, the area f Q h EBF (z;x) d e(z;x) 
is equal to zero. 

Theorem [2] combined with simulations (not included here) 
and the fact that an expected lower bound (analogous to 
Lemma 11 in | 6 ] ) on the number of independent guesses (in 
the asymptotic limit) for a Maxwell decoder computed from its 
DE recursion on the equivalent graph is equal to zero exactly 
at the area threshold (details omitted for brevity), leads to the 
following conjecture. 

Conjecture 1: The area threshold from Definition [2] e, is 
equal to the Maxwell decoding threshold and thus a lower 
bound on the MAP decoding threshold. 


We consider the ensemble (A ,p,N,w) of single-layer SC- 
CBs |j4l (coupling the original bipartite graph), where N is the 
coupling chain length and w, 1 < w < N + 1, is a smoothing 
parameter Q. The corresponding ensemble of SC graphs is 
denoted by Q C (A, p, N,w). The ensemble is constructed as 
follows. A collection of N flow-node groups are placed at 
positions {1,2 ,,N}, each containing k nodes of degree k, 
such that k = (the total number of flow nodes, i.e., the 
number of distinct flows to be counted, is Too). Furthermore, a 
collection of M = N + w — 1 counter-node groups are placed 
at positions {1,2 ,... ,M}, each containing - 

nodes of degree j (we implicitly assume that N is chosen such 
that 


K,Rn k 


and also '-jA are integers). 


The k k edge sockets in each group of flow and counter 
nodes are partitioned into w equally-sized subgroups (assum¬ 
ing that — is an integer) using a uniform random interleaver. 
We denote by and V^f\ the set of flow and counter node 
sockets, respectively, in the i-th subgroup, * = 0 , 1 ,..., w — 1 , 
at position n, where n = 1,... ,N for flow node sockets and 
n = 1 for counter node sockets. The SC ensemble 

is constructed by adding edges that connect the sockets in 
v\l\ to the sockets in i . Different graphs are obtained by 
different socket associations. Note that this construction leaves 
some sockets of the counter-node groups at the boundary 
unconnected and these are removed. In the following, we will 
denote the coupled ensemble Q C (A, p, N,w) by the alternative 
notation Q c (k,^,N,w). 


A. Density Evolution Recursion 

Denote by x) , i = 1,..., M , the error probability of an 
outgoing message from a flow node at position i at iteration 
£. Note that since there are no flow-node groups at positions 
i > N, initially x^ 1 = 0 for N < i < M. Using the ensemble 
defined above, we get the recursion shown in dl~ 2 l i at the top 
of the next page, where A = {A p q } is the N x M matrix 
defined by 


A 





if 1 < g — p+1 < w 
otherwise 


Note that contrary to SC-LDPC codes, for which the SC 
recursion contains two summations (an outer summation over 
N terms and an inner summation over M terms), the recursion 
for SC-CBs contains four summations (two summations over 
N terms and two summations over M terms), since the update 
rule for the flow nodes is different for odd and even iterations. 
Note that due to the band-structure of the matrix A there are 
only w nonzero terms in all four summations. Due to this 
difference between odd and even iterations, SC-CBs do not 
fit within the general framework of coupled scalar recursions 
outlined in 0 - , 

Let = (x{ ,..., x!"[l). The BP decoding threshold of 
the coupled ensemble Q c (k,"f,N,w) is defined (analogous to 
the uncoupled case) as e^p = d^ P (k, 7 , N, w), where 


eQ P (k,y, N,w) = sup je £ £ \ (k,"/,e, N,w) = oj . 
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Fig. 2. The difference e—£g P , where e is the area threshold from Definition^ 
and £g p is the BP decoding threshold of SC-CBs, as a function of fi for 
different values of the left-degree k with ( N,w) = (128,5) and (1,1) 
(uncoupled). 

VI. Numerical Results 

We have computed the BP decoding threshold ejjp for 
values of fi in the range 0.05 to 0.95 with a step size of 
0.05 for different values of the left-degree k for SC-CBs with 
(iV, w) = (128, 5) and (1,1) (uncoupled). For each value of /3 
(and k ) we have also computed the conjectured lower bound 
on the MAP decoding threshold (or the area threshold) from 
Conjecture[[] In Fig. [2] we plot the difference e e BP’ where e 
is the area threshold from Definition [2j as a function of /? = l ~. 
As we can observe from Fig. [2] there is a gap between the BP 
decoding threshold of SC-CBs (with ( N,w ) = (128,5)) and 
the conjectured lower bound on the MAP decoding threshold. 
However, the gap is significantly larger for uncoupled CBs, 
meaning that spatial coupling indeed improves performance. 
Note that the gap varies with fi and depends on k. 

VII. Discussion and Future Work 

As shown in Section [VI] we observe threshold improvement 
but not threshold saturation to the area threshold by coupling 
the original bipartite graph (see Section [VJ. We remark that 
coupling the equivalent graph representation mentioned in 
Section [III] gives exactly the same coupled DE recursion 
in (Ell. However, interestingly, if we consider the coupled 
version of ©-©, i.e., we substitute y and x in 0 and 0 , 
respectively, by an average over spatial positions, we do indeed 
observe threshold saturation to the area threshold. The DE 
recursion derived below corresponds to the coupled version of 
0 - 0 , split into two steps. 

Let a:) and y> , i = 1,, M, denote the output message 
error probability at a flow and counter node, respectively, at 
the £-th iteration at coupling chain position i. As before, we 


initialize xf = 0 for N < i < M. Now, define the following 
DE equations at iteration If, 

i-o(i-i eS? 1 -*'”-*’ x >-i n ) . is «*• 

1 — p ^1 — , if £ is even 

^ f (^'T 1 ■ ,sraM 

if « ^ even ' 

With this modified DE, threshold saturation to the area thresh¬ 
old which (from Conjecture 0 is a lower bound on the MAP 
decoding threshold, can be proved using the potential function 
framework by Yedla et al. outlined in 0. Note that this DE is 
characterized by an average only for odd and even iterations 
for the counter node and flow node updates, respectively. 
However, when coupling the original (or the equivalent) graph 
in the standard way, the average appears for all iterations 
(hence the four summations in (fTHl ). This effect, which is 
due to the fact that, as opposed to LDPC codes, the flow node 
update is different for odd and even iterations, seems to be the 
responsible for the lack of threshold saturation. 

A question that remains open is whether the DE equations 
above correspond to a physical system, i.e., whether threshold 
saturation can be achieved with an alternative coupling. Future 
research also includes a proof of Conjecture 1, the analysis 
of CBs with more layers, and exploring the connection with 
compressed sensing. 
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